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Abstract 

We study the dynamics of measure-valued solutions of what we call the 
EPDiff equations, standing for the Euler- Poincare equations associated with 
the diffeomorphism group (ofW 1 or an n-dimensional manifold M). Our main 
focus will be on the case of quadratic Lagrangians; that is, on geodesic motion 
on the diffeomorphism group with respect to the right invariant Sobolcv H 1 
metric. The corresponding Eulcr-Poincarc (EP) equations are the EPDiff equa- 
tions, which coincide with the averaged template matching equations (ATME) 
from computer vision and agree with the Camassa-Holm (CH) equations in one 
dimension. The corresponding equations for the volume preserving diffeomor- 
phism group are the well-known LAE (Lagrangian averaged Euler) equations 
for incompressible fluids. 

We first show that the EPDiff equations are generated by a smooth vector 
field on the diffeomorphism group for sufficiently smooth solutions. This is 
analogous to known results for incompressible fluids — both the Euler equations 
and the LAE equations — and it shows that for sufficiently smooth solutions, the 
equations are well-posed for short time. In fact, numerical evidence suggests 
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that, as time progresses, these smooth solutions break up into singular solutions 
which, at least in one dimension, exhibit soliton behavior. 

With regard to these non-smooth solutions, we study measure-valued so- 
lutions that generalize to higher dimensions the peakon solutions of the (CH) 
equation in one dimension. One of the main purposes of this paper is to show 
that many of the properties of these measure-valued solutions may be under- 
stood through the fact that their solution ansatz is a momentum map. Some 
additional geometry is also pointed out, for example, that this momentum map 
is one leg of a natural dual pair. 
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1 Introduction 

This paper is concerned with solutions of the EPDiff equations; that is, with the 
Euler-Poincare equations associated with the diffeomorphism group in n-dimensions. 
In particular, we are concerned with singular solutions that generalize the peakon 
solutions of the Camassa-Holm (CH) equation in one dimension. The CH equa- 
tion (see Camassa and Holm [1993]) for the dynamics of shallow water in a certain 
asymptotic regime, is 

u t + 3uu x = a 2 (u xxt + 2u x u xx + uu xxx ) , (1.1) 

where u(x, t) is the fluid velocity, subscripts denote partial derivatives in position x 
and time t, and a 2 is a positive constant. (See Dullin, Gottwald and Holm [2001, 
2003] for recent discussions of the derivation and asymptotic validity of the CH equa- 
tion for shallow water waves, at one order beyond the Korteweg-de Vries equation.) 
Equivalently, in Hamiltonian form, the CH equation reads 

rrit + um x + 2u x m = (1-2) 

where m = u — a 2 u xx and a 2 is a positive constant. As Camassa and Holm [1993] 
show, in Hamiltonian form the CH equation is the Lie-Poisson equation associated 
with the Lie algebra of one dimensional vector fields and with the Hamiltonian 

h(m) = — I umdx. (1-3) 
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The CH equation may be equivalently expressed in Euler-Poincare form by using 
the Lagrangian associated with the H 1 metric for the fluid velocity. That is, the 
Lagrangian as a function of the fluid velocity is given by the quadratic form, 

l(u) = \ f(u 2 + a 2 ul)dx. (1.4) 

It follows from Euler-Poincare theory (see Marsden and Ratiu [1999] and Holm, Marsden and Ratiu 
[1998a]) that the one parameter curve of diffeomorphisms r)(xo,t) depending on pa- 
rameter t and defined implicitly by 

d 

— 7](x ,t) = U(7](x ,t),t) 

is a geodesic in the group of diffeomorphisms of R (or, with periodic boundary 
conditions, of the circle S 1 ) equipped with the right invariant metric equal to the 
H 1 metric at the identity. 

A remarkable analytical property of the CH equation, conjectured by keep- 
ing track of derivative losses in Holm, Marsden and Ratiu [1998a] and proved in 
Shkoller [1998] is that the geodesic equations literally define a smooth vector field 
in the Sobolev H s topology. That is, in the material representation, the equa- 
tions have no derivative loss. This property is analogous to the corresponding 
results for the Euler equations for ideal incompressible fluid flow (discovered by 
Ebin and Marsden [1970]) and the Lagrangian averaged Euler equations (again con- 
jectured by Holm, Marsden and Ratiu [1998a] and proved by Shkoller [1998]). 

As we will explain in §3, a similar statement holds for the n-dimensional EPDiff 
equation if we use the H 1 metric. This is all the more remarkable because smooth- 
ness of the geodesic flow is not true for the 1? metric, at least not without assuming 
incompressibility. Smoothness of volume-preserving geodesic flow with respect to 
the I? metric does hold for the incompressible flow of an ideal Euler fluid, a result 
proved in Ebin and Marsden [1970]. 

Before proceeding with a discussion of the general case of the n-dimensional 
EPDiff equations, we shall quickly review, mostly to establish notation, a few facts 
about the Euler-Poincare and Lie-Poisson equations, whose basic theory is explained, 
for example, in Marsden and Ratiu [1999]. 

Review of Euler-Poincare and Lie-Poisson Equations. Let G be a Lie group 
and q its associated Lie algebra (identified with the tangent space to G at the identity 
element), with Lie bracket denoted by [£,rj\ for £, r/ € Q. Let £ : g — > R be a given 
Lagrangian. Let L : TG — > R be the right invariant Lagrangian on G obtained by 
translating £ from the identity element to other points of G via the right action 
of G on TG. A basic result of Euler-Poincare theory is that the Euler-Lagrange 
equations for L on G are equivalent to the (right) Euler-Poincare equations for £ on 
q, namely to 

d 51 ... 5£ 
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Here, adg : g — > g is the adjoint operator, the linear map given by the Lie bracket 
77 i— > [£,7)]. Also, ad^ : g* — > g* is its dual; that is, (ad|(/i), 77) = [£,77]), where 
( , ) is the natural pairing between g* and g. Also, 5£/5^ denotes the functional 
derivative of I with respect to £ £ g. For Ze/t invariant systems, we change the sign 
in (1.5). The Euler-Poincare equations can be written in the variational form 

!cft = 0, (1.6) 

for all variations of the form <5£ = 7) — [£, 77] for some curve 77 in g that vanishes at 
the endpoints. 

If the reduced Legendre transformation £ 1 — > /Li = 5£/5£ is invertible, then the 
Euler-Poincare equations are equivalent to the (right) Lie-Poisson equations: 

where the reduced Hamiltonian is given by, 

These equations are equivalent (via Lie-Poisson reduction and reconstruction) to 
Hamilton's equations on T*G relative to the Hamiltonian H : T*G — > M, obtained 
by right translating Zi from the identity element to other points via the right action 
of G on T*G. The Lie-Poisson equations may be written in the Poisson bracket 
form 

F = {F,h}, (1.8) 

where F : g* — > R is an arbitrary smooth function and the bracket is the (right) 
Lie-Poisson bracket given by 



{f,g}(m) = U 



5F_ SG 

S/j, ' S/j, 



(1.9) 



In the important case when £ is quadratic, the Lagrangian L is the quadratic 
form associated to a right invariant Riemannian metric on G. In this case, the Euler- 
Lagrange equations for L on G describe geodesic motion relative to this metric and 
these geodesies are then equivalently described by either the Euler-Poincare, or the 
Lie-Poisson equations. 



Outline of the paper. The main results of the present paper are as follows: 

1. In §2 we review some basic facts about the EPDiff equations, especially the 
singular solution ansatz (2.8) of Holm and Staley [2003] that introduces a class 
of singular solutions that generalize the peakon solutions of the CH equation 
to higher spatial dimensions. 

2. In §3 we give a plausibility argument that the EPDiff equations possess an in- 
teresting smoothness property; namely, they define a smooth vector field (that 
is, they define ODE's with no derivative loss) in the Lagrangian representation. 
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This means, in particular, that the EPDiff equations are locally well posed for 
sufficiently smooth initial data. Because of the development of singularities 
in finite time, which the numerics suggests, the smooth solutions may not 
exist globally in time. This smoothness property is similar to the correspond- 
ing smoothness property of the Euler equations for ideal incompressible fluid 
mechanics shown in Ebin and Marsden [1970]. 

3. In §4, we show that the singular solution ansatz (2.8) defines an equivariant 
momentum map. We do this in a natural way by identifying the singular 
solutions with certain curves in the space of embeddings Emb(<S, R n ) of a 
generally lower dimensional manifold S into the underlying manifold R n (or 
an n-manifold M) and letting the diffeomorphism group act on this space. 
The right action of Diff(S') corresponds to the right invariance of the EPDiff 
equations, while the left action of Diff(]R ri ) gives the desired solution ansatz. 

4. In §5 we explore the geometry of the singular solution momentum map, in 
parallel with the corresponding work on singular solutions (vortices, filaments, 
etc.) for the Euler equations of an ideal fluid in Marsden and Weinstein [1983]. 

5. Finally, in §6, we discuss some of the remaining challenges and speculate on 
some of the many possible future directions for this work. 

Historical Note. This paper is dedicated to our friend and collaborator Alan We- 
instein and, for us, this work parallels some of our earlier collaborations with him. 
Alan's basic works on reduction, Poisson geometry, semidirect product theory, and 
stability in mechanics — just to name a few areas — have been, and remain incredibly 
influential and important to the field of geometric mechanics. See, for instance, 
Marsden and Weinstein [1974]; Weinstein [1983b]; Marsden, Ratiu and Weinstein 
[1984]; Weinstein [1984]; Holm, Marsden, Ratiu and Weinstein [1985] 

Mechanics on Lie groups was pioneered by Arnold [1966], a reference that is a 
key foundation for the subject and in particular for the present paper. However, this 
theory was in a relatively primitive state, even by 1980, and it has benefited greatly 
from Alan's insights. In fact, the clear distinction between the Euler-Poincare and 
Lie- Poisson equations, the former with a variational structure and the latter with 
its Poisson structure, which took until the 1980's to crystalize, was greatly aided by 
Alan's work. 

Alan has made key contributions to many fundamental concepts in geometri- 
cal mechanics, such as Lagrangian submanifolds and related structures (Weinstein 
[1971, 1977]), symplectic reduction (Marsden and Weinstein [1974]), normal modes 
and periodic orbits (Weinstein [1973, 1978]), Poisson manifolds (Weinstein [1983b]), 
geometric phases (Weinstein [1990]), Dirac structures (Courant and Weinstein [1988]) 
groupoids and Lagrangian reduction (Weinstein [1996]) and the plethora of related 
"oid" structures he has been working on during the last decade (just look over the 
over 200 papers on MathSciNet he has written!) that will surely play an important 
role in the next generation of people working in the area of geometric mechanics. 
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Of Alan's papers, the one that is most directly relevant to the topics discussed in 
the present paper is Marsden and Weinstein [1983]. Alan himself is is still developing 
this area, too, as in Weinstein [2002]. 



2 The EPDiff Equation 

In this section we review the EPDiff equation; that is, the Euler-Poincare (EP) 
equation associated with the diffeomorphism group in ra-dimensions. This equa- 
tion coincides with a limiting case of the CH equation for shallow water waves in 
one and two dimensions. It also coincides with the ATME equation (the averaged 
template matching equation) in two dimensions. The latter equation arises in com- 
puter vision; see, for instance, Mumford [1998]; Hirani, Marsden and Arvo [2001] 
or Miller, Trouve and Younes [2002] for a description and further references. We 
have chosen to call this by a generic name, the EPDiff equation, because it has 
these various interpretations in different applications. Of course these different in- 
terpretations also provide opportunities: for example, to see to what extent the 
singular solutions found in the EPDiff equations are applicable, either for shallow 
water wave interactions, or for computer vision applications. A recent combination 
of these ideas in which image processing imitates soliton interactions appears in 
Holm, Trouve and Younes [2003]. 



Statement of the EPDiff Equations. Treating analytical issues formally at 
this point, let X denote the Lie algebra of vector fields on an n-dimensional man- 
ifold M (such as ]R n ). The vector fields comprise the algebra associated with the 
diffeomorphism group of M, but the usual Jacobi-Lie bracket is the negative of the 
(standard) Lie algebra bracket. (See Marsden and Ratiu [1999] for a discussion.) 

Let t : X — > R be a given Lagrangian and let 9Jt denote the space of one-form 
densities on M, that is, the momentum densities. The corresponding momentum 
density of the fluid is defined as 

m = — 6 9Jt, 
ou 

which is the functional derivative of the Lagrangian £ with respect to the fluid 
velocity If u is the basic dynamical variable, the EPDiff equations are simply 

the Euler-Poincare equations associated with this Lagrangian. Equivalently, if m 
is the basic dynamical variable, a Legendre transformation allows one to identify 
the EPDiff equations as the Lie-Poisson equations associated with the resulting 
Hamiltonian. For the case of R n , we will use vector notation for the momentum 
density m(x, t) : R n x R — > R n (a bold m instead of a lightface m). The EPDiff 
equations are as follows (see Holm, Marsden and Ratiu [1998a, b, 2002] for additional 
background and for techniques for computing the Euler-Poincare equations for field 
theories) , 

F) rp 

-m + u - Vrn + Vu ■ m + m (div u) = . (2-1) 



dt 

convection stretching expansion 
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In coordinates x % , i = 1,2, ... ,n, using the summation convention, and writing 
m = rriidx 1 (g> d n x (regarding m as a one-form density) and u = u l d/dx % (regarding 
u as a vector field), the EPDiff equations read 

d ■ drrii du J du J . N 

iJI'" ; • " 7hJ • m '7J7 • '"'^7 ,L (2 ' 2) 

The EPDiff equations can also be written very nicely as 

^ + £ u m = 0, (2.3) 

where £ u m denotes the Lie derivative of the momentum one form density m with 
respect to the velocity vector field u. 

As mentioned earlier, if I is a quadratic function of u, then the EPDiff equation 
((2.1) or, equivalently (2.3)), is the Eulerian description of geodesic motion on the 
diffeomorphism group of the underlying space (in this case M n ). The corresponding 
metric is the right invariant metric on the group whose value on the Lie algebra 
(the group's tangent space at the identity — the space of vector fields), is defined by 
I. Since the Lagrangian I is positive and quadratic in u, the momentum density is 
linear in u and so defines a positive symmetric operator Q op by 

m = — = <y op u . 

ou 

Variational Formulation. Following the variational formulation of EP theory, 
the particular EP equation (2.1) may be derived from the following constrained 
variational principle: 



5 

The variations are constrained to have the form 



; J £(u) dt = 0. 



5u = w + w • Vu — u • Vw. 



This assertion may of course be verified directly. These constraints are analogous 
to the so-called "Lin constraints" used for a similar variational principle for fluid 
mechanics. (See Marsden and Ratiu [1999] for a discussion and references.) 



Hamiltonian Formulation. A Legendre transformation yields the Hamiltonian, 

H{m) = (m,u) -£{u), 

where ( , ) is the natural pairing between one form densities and vector fields given 
by integration. This Hamiltonian is the corresponding quadratic form for the mo- 
mentum, 

ff(m)=*(Q-V)). (2.4) 
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Of course it often happens that Q op is a differential operator and in this case the 
inverse is usually given in terms of the convolution with the Green's function G, 
corresponding to the appropriate solution domain and boundary conditions; 

u = — = G * m . 

dm 

According to the general theory, the EP equation (2.1) may be expressed in 
Hamiltonian form by using the Lie-Poisson bracket on 9Jt as 

d 

— m = {m, H} LP = - ad * SH/Snx m. (2.5) 

One-dimensional CH Peakon Solutions. We return now to the CH equation 
(1.2), which, as we have noted, is the same as the EPDiff equation (2.1) for the case 
of one spatial dimension, when the momentum velocity relationship is defined by 
the Helmholtz equation, m = u — a 2 u xx . In one dimension, the CH equation has 
solutions whose momentum is supported at points on the real line via the following 
sum over Dirac delta measures, 

N 

m(x,t) = ^p i (t)6(x-q i (t)). (2.6) 

i=i 

The velocity corresponding to this measure-valued momentum is obtained by con- 
volution with the Green's function, 

G{\x-y\) = \e-\ x -y\' a : 

for the one-dimensional Helmholtz operator, Q op = (1 — a 2 d 2 ), appearing in the 
CH momentum velocity relationship, m = Q op u. Consequently, the CH velocity 
corresponding to this momentum is given by a superposition of peaked traveling 
wave pulses, 

1 N 

u{x,t) = -Y,Pi(t)^ qm/a - (2-7) 

i=l 

Thus, the superposition of "peakons" in velocity arises from the delta function 
solution ansatz (2.6) for the momentum. 

Remarkably, the isospectral eigenvalue problem for the CH equation implies 
that only these singular solutions emerge asymptotically in the solution of the ini- 
tial value problem in one dimension, Camassa and Holm [1993]. Figure 2.1 shows 
the emergence of peakons from an initially Gaussian velocity distribution and their 
subsequent elastic collisions in a periodic one-dimensional domain. 1 This figure 
demonstrates that singular solutions dominate the initial value problem and, thus, 
that it is imperative to go beyond smooth solutions for the CH equation; the situa- 
tion is similar for the EPDiff equation. 

1 Figure 2.1 was kindly supplied by Martin Staley 
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Figure 2.1: This figure shows a smooth localized (Gaussian) initial condition for the CH equation 
breaking up into an ordered train of peakons as time evolves (the time direction being vertical, 
which then eventually wrap around the periodic domain and interacting with other slower emergent 
peakons and causing a phaseshift (c.f. Alber and Marsden [1992]). 

Remarkably, the dynamical equations for pi(t) and qi(t), i = 1,...,N, that 
arise from solution ansatz (2.6-2.7) comprise an integrable system for any N. This 
system is studied in (Alber, Camassa, Fedorov, Holm and Marsden [2001]) and ref- 
erences therein. See also Vaninsky [2002, 2003] for discussions of how the integrable 
dynamical system for N peakons is related to the Toda chain with open ends. 

Generalizing the CH peakon solutions to n dimensions Building on the 
peakon solutions for the CH equation and the pulsons for its generalization to 
other traveling- wave shapes (see Fringer and Holm [2001]), Holm and Staley [2003] 
introduced the following measure-valued (that is, density valued) ansatz for the 
n— dimensional solutions of the EPDiff equation (2.1): 

N 

m(x, t)=J2 P a (s, t)S(x- Q a (s, t))ds. (2.8) 

a=l J 

These solutions are vector-valued functions supported in M n on a set of ./V surfaces 
(or curves) of codimension (n — k) for s€l' with k < n. They may, for example, 
be supported on sets of points (vector peakons, k = 0), one-dimensional filaments 
(strings, k = 1), or two-dimensional surfaces (sheets, k = 2) in three dimensions. 

One of the main results of this paper is that the singular solution ansatz (2.8) 
is a momentum map. This result helps to organize the theory and to suggest new 
avenues of exploration, as we shall explain. 

Substitution of the solution ansatz (2.8) into the EPDiff equations (2.1) implies 
the following integro-partial-differential equations (IPDEs) for the evolution of such 
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strings and sheets, 
8 N f 

¥t Q a (s,t) =J2 / P b (s',t)G(Q a (s,t)-Q b (s',t))ds', (2.9) 

6=1 J 

Ft N r 8 

^ a (s,t) = "Ey (r%s,tyP b (s\t))^^G(Q a (s,t)-Q\s',t))ds>. 

Importantly for the interpretation of these solutions, the coordinates s € R k turn out 
to be Lagrangian coordinates. The velocity field corresponding to the momentum 
solution ansatz (2.8) is given by 

N 

u(x,t) = G*m = Y^ / P b (s',t)G(x.-Q b {s',t))ds' , u£l". (2.10) 

6=1 

When evaluated along the curve x = Q a (s,t), the velocity satisfies, 
u(x,t)| / P b (s',t)G(Q a (s,t)-Q b ( s ',t))d S ' = ^^. (2.11) 

' 6=1 

Thus, the lower-dimensional support sets defined on x = Q a (s, t) and parameterized 
by coordinates s € R k move with the fluid velocity. Moreover, equations (2.9) for 
the evolution of these support sets are canonical Hamiltonian equations, 

1q-( s t ) - ^ —P a (s t)~- ^ (2 12) 

dt^ {S,t) ~ 5P" ' dt [S ' t} ~ SQ a ' [ ' 

The Hamiltonian function H N : (R™ x R n ) N -> R is, 

N 

Hn = II[J2 ( pa ( s ' *) • P '( s '' *)) G (Q a ( s > *) " Q 6 ( s '' *)) ds ds' . (2.13) 

^ a, 6=1 

This is the Hamiltonian for geodesic motion on the cotangent bundle of a set of 
curves Q a (s,t) with respect to the metric given by G. This dynamics was investi- 
gated numerically in Holm and Staley [2003] to which we refer for more details of 
the solution properties. 

One of our main goals is to show that the solution ansatz (2.8) can be phrased 
in terms of a momentum map that naturally arises in this problem. This geomet- 
ric feature underlies the remarkable reduction properties of the EPDiff equation 
and "explains" why the preceding equations must be Hamiltonian, namely, because 
momentum maps are Poisson maps. 

As explained in general terms in Marsden and Weinstein [1983], the way one 
implements a coadjoint orbit reduction is through a momentum map, and this holds 
even for the case of singular orbits (again ignoring functional analytic details). Thus, 
in summary, the reduction (2.12) is the EPDiff analog of the reduction in fluid 
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mechanics (that is, the EPDiffvoi equations) to point (or blob) vortex dynamics, 
vortex filaments, or sheets. 

There are, however, some important differences between vortex dynamics for in- 
compressible flows and the dynamics of the measure valued EPDiff solutions. For ex- 
ample, the Lagrangian representations of the equations of motion show that EPDiff 
solutions have inertia, while the corresponding solutions for point (or blob) vortices 
of the EPDiffvoi dynamics have no inertia. That is, the equations of motion for 
measure valued solutions on EPDiff yd are first order in time, while the dynamical 
equations for measure valued solutions on EPDiff are second order in time. This 
difference has profound effects on the properties of the solutions, especially on their 
stability properties. Numerical investigations of Holm and Staley [2003] show, for 
example, that the codimension-one solutions of EPDiff are stable, while higher codi- 
mension solutions of EPDiff are very unstable to codimension-one perturbations. In 
contrast, the co dimension- two solutions of EPDiffvd are known to be stable. 

Comments on the Physical Meaning of the Equations. The EPDiff equa- 
tions with the Helmholtz relation between velocity and momentum are not quite the 
CH equations for surface waves in 2D. Those would take precisely the same form, 
but the shallow water wave relation in the 2D CH approximation would be 

m = u — a 2 Grad Div u; that is, mi = U{ — a 2 ujji 
rather than the Helmholtz operator form, 

m = u — a 2 Div Gradti that is, mi = Ui — a 2 Uijj 
The corresponding Lagrangians are, respectively, 

hn(u) = ^J (\u\ 2 + a 2 (Divu) 2 )dxdy. (2.14) 

and 

Iepdis(u) = 2 J {\ u \ 2 + c* 2 \Gr&du\ 2 ) dx dy. (2.15) 

This difference was noted in Kruse, Scheurle and Du [2001], which identified 
(2.14) as the generalization of (1.4) for water waves in two dimensions. One may 
also verify this by considering the limit of the Green-Nagdhi equations for small 
potential energy. (The CH equation arises in this limit. The Lake and Great Lake 
equations of Camassa, Holm and Levermore [1996, 1997] also arise in a variant of 
this limit.) 

Remarkably, the numerics in Holm and Staley [2003] show that the solutions for 
a variety of initial conditions are indistinguishable in these two cases. The initial 
conditions in Holm and Staley [2003] were all spatially confined velocity distribu- 
tions. 

Notice that this difference affects the choice of Hamiltonian, but the equations 
are still Euler-Poincare equations for the diffeomorphism group and the description 
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of the anzatz (2.8) as a momentum map is independent of this difference in the 
equations. 

Figure 2.2 shows the striking reconnection phenomenon seen in the nonlinear 
interaction between wave-trains, as simulated by numerical solutions of the EPDiff 
equation and observed for internal waves in the Ocean. Fig 2.2(a) shows a frame 
taken from simulations of the initial value problem for the EPDiff equation in two di- 
mensions, excerpted from Holm and Staley [2003]. (See also Holm, Putkaradze and Stechmann 
[2003].) 

Fig 2.2(b) shows the interaction of two internal wave trains propagating at the 
interface of different density levels in the South China Sea, as observed from the 
Space Shuttle using synthetic aperture radar, courtesy of A. Liu (2002). Impor- 
tantly, both Fig 2.2(a) and Fig 2.2(b) show nonlinear reconnection occurring in the 
wave train interaction as their characteristic feature. 





Figure 2.2: Comparison of evolutionary EPDiff solutions in two dimensions (a) and Synthetic 
Aperture Radar observations by the Space Shuttle of internal waves in the South China Sea (b). 
Both Figures show nonlinear reconnection occurring in the wave train interaction as their charac- 
teristic feature. 

Internal waves and these sorts of interactions are generally thought to be de- 
scribed by the KP equation, and so any relations among the KP equation, the 
EPDiff equation and the 2D CH equation would be of great interest to explore; cf. 
Liu et al. [1998]. The derivations of the KP equation and the CH equation differ in 
the way the transverse motions are treated in the asymptotics; so some difference 
in their solution behavior is to be expected. 



3 Smoothness of the Lagrangian Equations 
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3 Smoothness of the Lagrangian Equations 

The One-Dimensional Case. Based on a formal argument given in Holm, Marsden and Ratiu 
[1998a], it was shown in Shkoller [1998] that when the CH equation (1.2) is trans- 
formed to Lagrangian variables, it defines a smooth vector field (that is, one obtains 
an evolution equation with no derivative loss). This means that one can show using 
ODE methods that the initial value problem is well-posed and one may also establish 
other important properties of the equations when the data is sufficiently smooth. 

As above, we write the relation between m and u as m = Q op u, so that in the 
one dimensional case, Q op is the operator Q op = Id — a 2 d xx . We first recall how the 
equations are transformed into Lagrangian variables. Introduce the one parameter 
curve of diffeomorphisms rf(xo,t) defined implicitly by 

d 

— rj(x , t) = u(rj(x ,t),t), (3.1) 

so that r] is a geodesic in the group of diffeomorphisms of R (or, with periodic 
boundary conditions, of the circle S 1 ) equipped with the right invariant metric 
equaling the H 1 metric at the identity. 

We compute the second time derivative of rj in a straightforward way by differ- 
entiating (3.1) using the chain rule: 



d 2 r\ du 
-— - = uu r + — — 
dt 2 dt 



Acting on this equation with Q op and using the definition m = Q op u yields 

d 2r rj din 
Qop -q£2 = Q°p( uu z) - ud x (Q op u) + um x + — 

= [Qo P , ud x ]u - 2mu x 
= —3a u x u xx — 2mu x 

where the third step uses the commutator relation calculated from the product rule, 

[Q op , ud x ]u = -3a 2 u x u xx . 
Hence, the preceding equation becomes 

d 2 r] 1 



dt 2 ~ 2 



Qop 1 dx [p?u\ + 2u 2 ) . (3.2) 



The important point about this equation is that the right hand side has no 
derivative loss. That is, if u is in the Sobolev space H s for s > 5/2, then the right 
hand is also in the same space. Regarding the right hand side as a function of rj 
and drj/dt, we see that it is plausible that the second order evolution equation (3.2) 
for rj defines a smooth ODE on the group of H s diffeomorphisms. (This argument 
requires the use of, for example, weighted Sobolev spaces in the case x € R). 
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The above is the essence of the argument given in Shkoller [1998], Remark 3.5, 
which in turn makes use of the type of arguments found in Ebin and Marsden [1970] 
for the incompressible case and which shows, by a more careful argument, that the 
spray is smooth if s > 3/2. However, one should note that the complete argument 
is not quite so simple (just as in the case of incompressible fluids). A subtilty arises 
because smoothness means as a function of rj , r). Hence, one must express u in terms 
of rj, namely through the relation m = rjt ° Vt^ 1 > where the subscript t here denotes 
that this argument is held fixed, and is not a partial derivative. Doing this, one sees 
that, while there is clearly no derivative loss, the right hand side of (3.2) does involve 
r^ 1 and the map rjt >— > is known to not be smooth (just as in Ebin and Marsden 
[1970]). Nevertheless, the combination that appears in (3.2) is, quite remarkably, a 
smooth function of 77, r). 



The ra-Dimensional Case. The above argument readily generalizes to the n- 
dimensions, which we shall present in the case of W 1 or the flat n-torus T n for 
simplicity. Namely, we still have the relation 

d 

— ?7(xo,t) = u(r?(x ,t),t), 

between r\ and u. Consequently, we may compute the second partial time derivative 
of rj in the usual fashion using the chain rule: 

d 2 r\ du 

Therefore, as in the one dimensional case, we get 

Qop-rp" = [Qop , (u • V)] u + u • Vm + — . 
Calculating the commutator relation in n-dimensions gives 

[Q op , (u • V)] u = - a 2 div ( Vu • Vu + Vu • Vu T ) + a 2 ( Vu) • Vdiv u 
or, in components, 

( [Q op , (u • V)]u) i = - a 2 d k (u it j u jtk + u itj u k ,j) + a 2 (u itj ) djdivu , (3.3) 

with a sum on repeated indices. 

Upon substituting the preceding commutator relation, the EPDiff equation (2.1) 
and the vector calculus identity 

-Vu r • m = a 2 div (Vu T • Vu) - V^|u| 2 + y|Vu| 2 ) , (3.4) 
then imply the n— dimensional result 

^ ^1 „,2j;, r fv7,, V7,, 1 V7,, V7,,T V7.,T V7,, V7,,7Vj- „\ , ^T^|V7,,|2 



Qop-gjj = -a 2 div(^Vu • Vu + Vu • Vu J - Vu 7 • Vu - Vu 1 (div u) + - Id|Vu 
— u(divu) — - V|u| 2 . 
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This form of the EPDiff equation is useful for interpreting some of its solution 
behavior. As in the one dimensional case, the crucial point is that the right hand 
side involves at most second derivatives of u; so there is no derivative loss in the 
overall expression. The precise formulation of the smoothness result thus should 
hold, as in the one-dimensional case, although, naturally the details are a bit more 
complicated. We shall leave the technical details for another publication. 

A consequence is that the equations are well posed for short time if the initial 
data is smooth enough and, hence, e.g., two nearby smooth solutions can be joined by 
a unique geodesic. Moreover, the solutions of the EPDiff equation are automatically 
C°° in time. In the sections that follow, we will be interested in nonsmooth data, 
which is in stark contrast to the preceding discussion, which requires initial data 
that is at least C . 

Remarkably, the same smoothness results hold for the case of the LAE-a (La- 
grangian averaged Euler) equations, a set of incompressible equations in which 
small scale fluctuations are averaged. One can view the LAE-a equations as the 
incompressible version of the EPDiff equations. This smoothness property for 
the LAE-a equations was shown by Shkoller [1998] for regions with no boundary 
and for regions with boundary (for various boundary conditions), it was shown in 
Marsden, Ratiu and Shkoller [2000]. However, unlike the incompressible case, the 
results apparently do not hold if a is zero (as also noted in Shkoller [1998]). This 
sort of smoothness result also appears not to hold for many other equations, such as 
the KdV equation, even though it too can be realized as Euler-Poincare equations 
on a Lie algebra, or as geodesies on a group, in this case the Bott-Virasoro group, 
as explained in Marsden and Ratiu [1999] and references therein. 

The Development of Singularities. The smoothness property just discussed 
does not preclude the development in finite time of singular solutions from smooth 
localized initial data as was indicated in Figure 2.1. To capture the local singularities 
in the EPDiff solution (either verticality in slope, or discontinuities in its spatial 
derivative) that develop in finite time from arbitrarily smooth initial conditions, 
one must enlarge the solution class of interest, by considering weak solutions. 

There are a number of papers on weak solutions of the CH equation such as 
Xin and Zhang [2000] that we will not survey here. We just mention that the 
theory is not yet complete, as it is still unknown in what sense one may de- 
fine global unique weak solutions to the CH equations in H . As discussed in 
Alber, Camassa, Fedorov, Holm and Marsden [2001] for the CH equation, one most 
likely must consider weak solutions in the spacetime sense. 

The steepening lemma of Camassa and Holm [1993] proves that in one dimen- 
sion any initial velocity distribution whose spatial profile has an inflection point 
with negative slope (for example, any antisymmetric smooth initial distribution of 
velocity on the real line) will develop a vertical slope in finite time. Note that the 
peakon solution (2.7) has no inflection points, so it is not subject to the steepening 
lemma. However, the steepening lemma underlies the mechanism for forming these 
singular solutions, which are continuous but have discontinuous spatial derivatives; 
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they also lie in H l and have finite energy. We conclude that solutions with initial 
conditions in H s with s > (n/2) + 1 go to infinity in the H s norm in finite time, 
but remain in H 1 and presumably continue to exist in a weak spacetime sense for 
all time in H 1 . 

Numerical evidence in higher dimensions and the inverse scattering solution for 
the CH equation in one dimension (the latter has only discrete eigenvalues, corre- 
sponding to peakons) both show that the singular solutions completely dominate 
the time-asymptotic dynamics of the initial value problem (IVP). This singular IVP 
behavior is one of the main discoveries of Camassa and Holm [1993]. This singular 
behavior has drawn a great deal of mathematical interest to the CH equation and 
its relatives, such as EPDiff. The other properties of CH — its complete integra- 
bility, inverse scattering transform, connections to algebraic geometry and elliptical 
billiards, bi-Hamiltonian structure, etc. - are all interesting, too. However, the 
requirement of dealing with singularity as its main solution phenomenon is the pri- 
mary aspect of CH (and EPDiff). We aim to show that many of the properties 
of these singular solutions of CH and EPDiff are captured by recognizing that the 
singular solution ansatz itself is a momentum map. This momentum map property 
explains, for example, why the singular solutions (2.8) form an invariant manifold 
for any value of N. 

In one dimension, the complete integrability of the CH equation as a Hamilto- 
nian system and its soliton paradigm completely explain the emergence of peakons 
in the CH dynamics. Namely, their emergence reveals the initial condition's soli- 
ton (peakon) content. However, beyond one dimension, we do not have an explicit 
mechanism for explaining why only singular solution behavior emerges in numerical 
simulations. One hopes that eventually a theory will be developed for explaining 
this singular solution phenomenon in higher dimensions. Such a theory might, for 
example, parallel the well-known explanation of the formation of shocks for hyper- 
bolic partial differential equations. (Note, however, that EPDiff is not hyperbolic, 
because the relation u = G * m between its velocity and momentum is nonlocal.) 

In the remainder of this work, we shall focus our attention on the momentum 
map properties of the invariant manifold of singular solutions (2.8) of the EPDiff 
equation. 

4 The Singular Solution Momentum Map 

The Momentum Ansatz (2.8) is a Momentum Map. The purpose of this 
section is to show that the solution ansatz (2.8) for the momentum vector in the 
EPDiff equation (2.1) defines a momentum map for the action of the group of 
diffeomorphisms on the support sets of the Dirac delta functions. These support 
sets are the analogs of points on the real line for the CH equation in one dimension. 
They are points, curves, or surfaces in M n for the EPDiff equation in n— dimensions. 

This result, as we shall discuss in greater detail later, shows that the solution 
ansatz (2.8) fits naturally into the scheme of Clebsch, or canonical variables in the 
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sense advocated by Marsden and Weinstein [1983] as well as showing that these 
singular solutions evolve on special coadjoint orbits for the diffeomorphism group. 

One can summarize by saying that the map that implements the canonical (Q, P) 
variables in terms of singular solutions is a (cotangent bundle) momentum map. 
Such momentum maps are Poisson maps; so the canonical Hamiltonian nature of 
the dynamical equations for (Q, P) fits into a general theory which also provides a 
framework for suggesting other avenues of investigation. 

Theorem 4.1. The momentum ansatz (2.8) for measure-valued solutions of the 
EPDiff equation (2.1), defines an equivariant momentum map 

Jsing : r*Emb(S,R n ) -» X(R n )* 

that we will call the singular solution momentum map. 

We shall explain the notation used in this statement in the course of the proof. 
Right away, however, we note that the sense of "defines" is that expressing m in 
terms of Q, P (which are, in turn, functions of s) can be regarded as a map from the 
space of (Q(s), P(s)) to the space of m's. This will turn out to be the Lagrange-to- 
Euler map for the fluid description of the singular solutions. 

We shall give two proofs of this result from two rather different points of view. 
The first proof below uses the formula for a momentum map for a cotangent lifted 
action, while the second proof focuses on a Poisson bracket computation. Each proof 
also explains the context in which one has a momentum map. (See Marsden and Ratiu 
[1999] for general background on momentum maps.) 

First Proof. For simplicity and without loss of generality, let us take N = 1 and so 
suppress the index a. That is, we shall take the case of an isolated singular solution. 
As the proof will show, this is not a real restriction. 

To set the notation, fix a fc-dimensional manifold S with a given volume element 
and whose points are denoted s € S. Let Emb(5, M n ) denote the set of smooth em- 
beddings Q : S — » R n . (If the EPDiff equations are taken on a manifold M, replace 
M n with M.) Under appropriate technical conditions, which we shall just treat for- 
mally here, Emb(S l , M n ) is a smooth manifold. (See, for example, Ebin and Marsden 
[1970] and Marsden and Hughes [1983] for a discussion and references.) 

The tangent space Tq Emb(5, IR n ) to Emb(S,R n ) at the point Q € Emb(5',R n ) 
is given by the space of material velocity fields, namely the linear space of maps 
V : S — > W 1 that are vector fields over the map Q. The dual space to this space will 
be identified with the space of one- form densities over Q, which we shall regard as 
maps P : S — > (R n )*. In summary, the cotangent bundle T* Emb(S f , MJ 1 ) is identified 
with the space of pairs of maps (Q, P). 

These give us the domain space for the singular solution momentum map. Now 
we consider the action of the symmetry group. Consider the group (5 = Diff of 
diffeomorphisms of the space S in which the EPDiff equations are operating, con- 
cretely in our case M. n . Let it act on © by composition on the left. Namely for 
n G Diff (IT), we let 

7/-Q = t/oQ. (4.1) 
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Now lift this action to the cotangent bundle T* Emb(5, W 1 ) in the standard way (see, 
for instance, Marsden and Ratiu [1999] for this construction). This lifted action is 
a symplectic (and hence Poisson) action and has an equivariant momentum map. 
We claim that this momentum map is precisely given by the ansatz (2.8). 

To see this, we just need to recall and then apply the general formula for the mo- 
mentum map associated with an action of a general Lie group (5 on a configuration 
manifold Q and cotangent lifted to T*Q. 

First let us recall the general formula. Namely, the momentum map is the map 
J : T*Q — > g* (q* denotes the dual of the Lie algebra g of (5) defined by 

J(a,)-£ = (a g ,fo(g)>, (4.2) 

where a q £ T*Q and £ 6 Q, where £q is the infinitesimal generator of the action of 
on Q associated to the Lie algebra element £, and where (a q ,^Q(q)) is the natural 
pairing of an element of T*Q with an element of T q Q. 

Now we apply this formula to the special case in which the group © is the 
diffeomorphism group Diff(IR n ), the manifold Q is Emb(5, R n ) and where the action 
of the group on Emb(5, M n ) is given by (4.1). The sense in which the Lie algebra 
of © = Diff is the space q = X of vector fields is well-understood. Hence, its dual is 
naturally regarded as the space of one- form densities. The momentum map is thus 
a map J : T* Emb(5,R n ) -► X*. 

With J given by (4.2), we just need to work out this formula. First, we shall 
work out the infinitesimal generators. Let X £ X be a Lie algebra element. By 
differentiating the action (4.1) with respect to 77 in the direction of X at the identity 
element we find that the infinitesimal generator is given by 

^Emb(5,R n )(Q) = X o Q. 

Thus, taking ot q to be the cotangent vector (Q,P), equation (4.2) gives 

(J(Q,P),X) = ((Q,P),XoQ) 

= f P 4 ( S )X 4 (Q( s ))d fc s. 
Js 

On the other hand, note that the right hand side of (2.8) (again with the index a 
suppressed, and with t suppressed as well), when paired with the Lie algebra element 
X is 

l^j P(s)5(x-Q(s))d k s,X^ = J J (Pi(s)S(x- Q(s))d k s^ X l (x)d n x 

= / P i {s)X i (Q(s)d k s. 
Js 

This shows that the expression given by (2.8) is equal to J and so the result is 
proved. ■ 
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Second Proof. As is standard (see, for example, Marsden and Ratiu [1999]), one 
can characterize momentum maps by means of the following relation, required to 
hold for all functions F on T* Emb(»S, ]R n ); that is, functions of Q and P: 

{F,{J,Z)} = Zp[F]. (4.3) 

In our case, we shall take J to be given by the solution ansatz and verify that it 
satisfies this relation. To do so, let £ G X so that the left side of (4.3) becomes 

On the other hand, one can directly compute from the definitions that the infinitesi- 
mal generator of the action on the space T* Emb(S l , M. n ) corresponding to the vector 
field £*(x)g|p- (a Lie algebra element), is given by (see Marsden and Ratiu [1999], 
formula (12.1.14)): 

which verifies that (4.3) holds. 

An important element left out in this proof so far is that it does not make clear 
that the momentum map is equivariant, a condition needed for the momentum map 
to be Poisson. The first proof took care of this automatically since momentum maps 
for cotangent lifted actions are always equivariant and hence Poisson. 

Thus, to complete the second proof, we need to check directly that the mo- 
mentum map is equivariant. Actually, we shall only check that it is infinitesimally 
invariant by showing that it is a Poisson map from T*Emb(5, W 1 ) to the space of 
m's (the dual of the Lie algebra of 3C) with its Lie- Poisson bracket. This sort of 
approach to characterize equivariant momentum maps is discussed in an interesting 
way in Weinstein [2002]. 

The following computation accomplishes this methodology by showing that the 
singular solution momentum map is Poisson. 

Indeed, we use the canonical Poisson brackets for {P}, {Q} and apply the chain 
rule to compute {mj(x), mj(y)}, with notation 5' k (y) = dS(y)/dy k . We get 



S nQ(5)) 



5F 5 
5Pj 5Qi 



■ems)) 



d k s 



(m i (x),m i (y)} 



N N . 

ds p a (s, t) «J(x - Q a (s, t)) , Y, ds'P^s 1 , t) S(y - Q b (s', t)) 
a =r b=l J J 



N 



,,6=1 



^2 / jdsds' 



{P?(s), P^s')} *(x - Q a (s)) f (y - Q V)) 



bfj\ 



- {P?(s), Q b k (s')}P>(s') «5(x - Q a (s)) 5 > k (y - Q b (s')) 

- {Qt(s),P!(s')}P?(s)6' k (x - Q%s)) «5(y - Q b (s')) 

+ {Qt(s), Q b e (s')}P t a (s)PHs')5' k (x ~ Q a (s)) 5'tiy - Q b (s')) 
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Substituting the canonical Poisson bracket relations 

{jf( s ),pjV)} = o 

{Q a k (s),Q b e (s')} = 0, and 
{QUs),P^s')} = 5 ab 5 kj 5(s 

into the preceding computation yields, 

{m i (x),m i (y)} 

N „ TV 



a=:r 6=1^ 
TV 

£ /d S ^(,)5(x-Q a ( S ))5:(y-Q a ( S )) 

a=l J 



N 

/^( S )5 / 7 (x-Q a ( S ))5(y-Q a ( S )) 

o=l •* 



3 

d 



^ (x te + ^7 mi(x) j 5(x - y) - 



Thus, 



which is readily checked to be the Lie-Poisson bracket on the space of m's. This 
completes the second proof of theorem. ■ 

Each of these proofs has shown the following basic fact. 

Corollary 4.2. The singular solution momentum map defined by the singular so- 
lution ansatz, namely, 

Jsing : T* Emb(S,R n ) -» X(R n )* 

is a Poisson map from the canonical Poisson structure on T* Emb(5, W 1 ) to the 
Lie-Poisson structure on X(M n )* . 

This is perhaps the most basic property of the singular solution momentum map. 
Some of its more sophisticated properties are outlined in the following section. 



Pulling Back the Equations. Since the solution ansatz (2.8) has been shown in 
the preceding Corollary to be a Poisson map, the pull back of the Hamiltonian from 
X* to T*Emb(5, W 1 ) gives equations of motion on the latter space that project to 
the equations on X*. This is why the functions Q a (s, t) and P a (s, t) satisfy canonical 
Hamiltonian equations. Note that the coordinate s G M fc that labels these functions 
is a Lagrangian coordinate. 
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In terms of the pairing 

(-,■}: *X0^R, (4.5) 

between the Lie algebra q (vector fields in M. n ) and its dual q* (one- form densities in 
W 1 ), the following relation holds for measure- valued solutions under the momentum 
map (2.8), 

(m , u) = J m • ud n x, L 2 pairing for m&u £ ]R n , 
N 



II {P a (s,t)-P b (s',t))G{Q a (s,t)-Q b (s',t))dsds' 
dQ a {s,t) 



a,b=l 
N 

/^P a (s,t) • 



dt 

a=l 



= «P,Q», (4-6) 

which is the natural pairing between the points (Q, P) € T* Emb(5, M n ) and (Q, Q) € 
TEmb(S, R n ). 

The pull-back of the Hamiltonian i?[m] defined on the dual of the Lie algebra 
q*, to T* Emb(S', M n ) is easily seen to be consistent with what we had before: 

H[m] = i(m , G * m) = i«P , G * P» = H N [P, Q] . (4.7) 

In summary, in concert with the Poisson nature of the singular solution momen- 
tum map, we see that the singular solutions in terms of Q and P satisfy Hamiltonian 
equations and also define an invariant solution set for the EPDiff equations. In fact, 
this invariant solution set is a special coadjoint orbit for the diffeomorphism group, 
as we shall discuss in the next section. 



Smoothness. It would be extremely interesting if the smoothness properties ex- 
plored in §3 were also valid on the space T* Emb(S l , W 1 ). We hope to explore this 
point in future publications. 



5 The Geometry of the Momentum Map 

In this section we explore the geometry of the singular solution momentum map 
discussed in §4 in a little more detail. The main idea may be stated as follows: 
simply apply all of the ideas given in Marsden and Weinstein [1983] in a systematic 
way to the current setting. 



Coadjoint Orbits. The first claim is that the image of the singular solution mo- 
mentum map is a coadjoint orbit in X* . This means that (perhaps modulo some 
issues of connectedness and smoothness, which we do not consider here) the solution 
anzatz given by (2.8) defines a coadjoint orbit in the space of all one-form densities, 
regarded as the dual of the Lie algebra of the diffeomorphism group. 
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These coadjoint orbits should be thought of as singular orbits — that is, due to 
their special nature, they are not generic. Also, this explains why the singular 
solutions form dynamically invariant sets — it is because they are coadjoint orbits, 
which are symplectic submanifolds of the Lie-Poisson manifold 3£(R n )*. 

The idea of the proof of this is simply this: whenever one has an equivariant 
momentum map J : P — > 0* for the action of a group G on a symplectic or Poisson 
manifold P, and that action is transitive, then the image of J is an orbit (or at least 
a piece of an orbit). This general result, due to Kostant, is stated more precisely 
in Marsden and Ratiu [1999], Theorem 14.4.5. Roughly speaking, the reason that 
transitivity holds in our case is because one can "move the manifolds S around at 
will" using diffeomorphisms. 

Symplectic Structure on Orbits. Recall (from, for example, Marsden and Ratiu 
[1999]), the general formula for the symplectic structure on coadjoint orbits: 



where fx € Q* is a chosen point on an orbit and where £, rj are elements of g. We use 
a plus sign in this formula since we are dealing with orbits for the right action. 

Just as in Marsden and Weinstein [1983], this leads to an explicit formula for 
the coadjoint orbit symplectic structure in the case of Diff . In the present case, it 
is a particularly simple and transparent formula. 

Recall that in the case of incompressible fluid mechanics, this procedure leads 
naturally to the symplectic (and Poisson) structure for many interesting singular 
coadjoint orbits, such as point vortices in the plane, vortex patches, vortex blobs 
(closely related to the planar LAE-a equations) and vortex filaments. 

For the case of the diffeomorphism group, let O m denote the coadjoint orbit 
through the point m G JC*(R n ). 

Theorem 5.1. The symplectic structure £l m on T m O m is given by 



Proof. We simply plug into the general Kirillov-Kostant-Souriau formula (5.1) for 
the symplectic structure on coadjoint orbits. (As noted above, there is a + sign, since 
we are dealing with a right invariant system) . The only thing needing explanation is 
that our Lie algebra convention always uses the left Lie bracket. For Diff, this is the 
negative of the usual Lie bracket as is explained in Marsden and Ratiu [1999]. ■ 

Dual Pairs. The singular solution momentum map Jsing : T* Emb(5, W 1 ) — > 
X(W n )* forms one leg of a dual pair. The point is that there is another group 
that acts on Emb(5, W n ), namely the group Diff (5) of diffeomorphisms of S, which 
acts on the right, while Diff (BL n ) acted on it by composition on the left (and which 
gave rise to our singular solution momentum map). This momentum map will be 



(£b*(m),%-(^)) = <A*> [£>»7]> 



(5.1) 
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denoted Js '■ T* Emb(5, M n ) — > 3C(S)* . The dual pair arises for reasons that are 
very similar to those in Marsden and Weinstein [1983]. See the following diagram. 



We also note that this framework shows that the parameterization of the singular 
solutions in terms of Q and P are exactly Clebsch variables in the sense given in 
Marsden and Weinstein [1983]. Also notice that when we write the singular solutions 
in Q-P space, we are finding solutions that are collective and so all the properties 
of collectivization are valid. See Marsden and Ratiu [1999] for a general discussion 
and references to the original work of Guillemin and Sternberg on this topic. 

6 Challenges, Future Directions and Speculations 

Numerical Issues: Geometric Integrators. The computations of Martin Stan- 
ley that were illustrated in this paper and that are discussed in Holm and Staley 
[2003], make use of both mimetic differencing and reversibility in a critical way 
and this is important for good numerical resolution. In other words, integra- 
tors that respect the basic geometry underlying the problem seem to play an key 
role. It would be interesting to pursue this aspect further and also incorporate 
discrete exterior calculus and variational multisymplectic integration methods (see 
Desbrun, Hirani, Leok and Marsden [2003] as well as Marsden, Patrick and Shkoller 
[1998] and Lew, Marsden, Ortiz and West [2003]). 

Analytical Issues: Geodesic Incompleteness of H 1 EPDiff. The emergence 
in finite time of singular solutions from smooth initial data observed numerically 
in Holm and Staley [2003] indicates that the diffeomorphism group with respect to 
the right invariant H 1 metric is geodesically incomplete when the diffeomorphism 
group has the H s topology, s > (n/2) + 1. The degree of its geodesic incompleteness 
is not known, but we suspect that almost all EPDiff geodesies in H 1 cannot be 
extended indefinitely. This certainly holds in one spatial dimension, where the 
discreteness of the CH isospectrum implies that asymptotically in time the CH 
solution arising from any confined initial velocity data consists only of peakons. It 
is an important challange to find a context in which one can put the H 1 topology 
on the diffeomorphism group and reestablish geodesic completeness. The numerics 
suggests that this might be possible, while known existence theorems, even for the 
CH equation are not yet capable of showing this — to the best of our knowledge. 



T* Emb(5, W 1 ) 




x(R n y 



x(sy 
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Reversible Reconnections of the Singular EPDifF Solutions. EPDiff is a 
reversible equation, and the collisions of its peakon solutions on the line 1R 1 (or the 
circle S 1 ) are known to be reversible. In principle, the reconnections of the sin- 
gular EPDiff solutions observed numerically in Holm and Staley [2003] in periodic 
domains T 2 and T 3 must also be reversible. Reversibility of its reconnections dis- 
tinguishes the singular solutions of EPDiff from vortex fluid solutions and shocks in 
fluids, whose reconnections apparently require dissipation and so, are not reversible. 
The mimetic finite differencing scheme used for the numerical computation of EPDiff 
solutions in Holm and Staley [2003] was indeed found to be reversible for overtaking 
collisions, but it was found to be only approximately reversible for head-on collisions, 
which are much more challenging for numerical integration schemes. 

Applications of EPDiff Singular Solutions in Image Processing. The sin- 
gular EPDiff solutions correspond to outlines (or cartoons) of images in applications 
of geodesic flow for the template, or pattern matching approach. The dynamics of 
the singular EPDiff solutions described by the momentum map (2.8) introduces the 
paradigm of soliton collisions into the mechanics and analysis of image processing 
by template matching. (See Holm, Trouve and Younes [2003] for more discussions 
of this new paradigm for image processing.) The reversibility of the collisions among 
singular solutions and their reconnections under EPDiff flow assures the preserva- 
tion of the information contained in the image outlines. In addition, the invariance 
of the manifold of iV singular solutions under EPDiff assures that the fidelity of 
the image is preserved in the sense of approximation theory. That is, an iV soliton 
approximation of the image outlines remains so, throughout the EPDiff flow. A 
natural approach for numerically simulating EPDiff flows in image processing is to 
use multisymplectic algorithms. The preservation of the space-time multisymplectic 
form by these algorithms introduces an initial- value, final-value formulation of the 
numerical solution procedure that is natural for template matching. 

Rigorous Poisson Structures. In Vasylkevych and Marsden [2003], the ques- 
tion of the (rigorous) Poisson nature of the time t map of the flow of the Euler 
equations for an ideal fluid in appropriate Sobolev spaces is explored. Given the 
smoothness properties in §3, it seems reasonable that similar properties should also 
hold for the EPDiff equations. However, as mentioned earlier, these smoothness 
properties do not preclude the emergence of singular solutions from smooth initial 
data in finite time, because of the possibility for geodesic incompleteness. 

Other Groups. The general setting of this paper suggests that perhaps one should 
look for similar measure valued or singular solutions associated with other problems, 
including geodesic flows on the group of symplectic diffeomorphisms (relevant for 
plasma physics, as in Marsden and Weinstein [1982]), Bott-Virasoro central exten- 
sions and super-symmetry groups. 

Scattering. It might be interesting to explore the relation of the singular solu- 
tion momentum map (2.8) to integrability and scattering data. For example, see 
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Vaninsky [2003] for an interesting discussion of the Poisson bracket for the scatter- 
ing data of CH in ID. This turns out to be the Atiyah-Hitchin bracket, which is 
also related to the Toda lattice, and this fascinating observation leads to an infinite- 
dimensional version of Jacobi elliptic coordinates. 

Other Issues. Of course there are many other issues remaining to explore that 
are suggested by the above setting, such as convexity of the momentum map, its 
extension to Riemannian manifolds, etc. We shall, however, leave these issues for 
other publications and other researchers. 
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